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We present a detailed analytical and numerical analysis of the nuclear spin dynamics in parabolic 
quantum wells. The shallow potential of parabolic quantum wells permits substantial modification 
of the electronic wave function in small electric fields. The nuclear spin relaxation via the hyperfine 
interaction depends on the electronic local density of states, therefore the local nuclear relaxation 
time depends sensitively on the electric field. For an inhomogeneous nuclear magnetization, such as 
generated by dynamic nuclear polarization, the total nuclear magnetization dynamics can similarly 
be altered. We examine this effect quantitatively and the effect of temperature, field, well thickness, 
and nuclear spin diffusion. 
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I. INTRODUCTION 

The long electronic spin coherence times in GaAs 
structures [1] may assist the practical realization of 
spin based electronics, including quantum information 
processing. [2] Electronic spin interactions with phonons, 
nuclear spins, other electron spins and impurities will 
limit the value of the electronic spin relaxation and coher- 
ence times. Due to much longer coherence times nuclear 
spins can also be useful elements of spin-based devices. 
Compared to electronic spin relaxation times, of the or- 
der of 100 ns, nuclear spin relaxation times in the absence 
of carriers are roughly 10 min in GaAs structures. [3-5] 
Yet the local electronic density of states near a nuclear 
spin can be increased to the point where the hyperfine 
interaction dominates nuclear spin relaxation and the re- 
laxation times are shorter than 1 s. 

In this work we investigate the nuclear spin relaxation 
time in a parabolic quantum well (PQW). In this low di- 
mensional structure electrons are trapped in a parabolic 
potential that restricts their motion along one direction. 
Such a profile for the confining potential can be ob- 
tained by gradually varying the Al concentration, x, of 
Al^Gai^^As layers between two barriers. [6] In the pres- 
ence of an electric field the minimum of the confining 
potential is displaced with respect to the center of the 
PQW, whereas the shape of the confining potential and 
implicitly of the envelope wave function are largely con- 
served. This particular property of the PQW makes it 
suitable for the practical realization of electronic devices 
based on the electrical manipulation of the electronic or 
nuclear spins. Recent experiments [7] demonstrated that 
considerable electrical control of electronic spin dynamics 
(g-factors) can be achieved in shallow potential parabolic 
quantum wells. 

The standard technique used to extract nuclear spin 
relaxation times is nuclear magnetic resonance (NMR) 
spectroscopy. For low dimensional systems such as 
PQW's the signal in a standard NMR experiment from 
the number of available nuclei across the sample is too 
small to measure. The acquired polarization of nu- 
clei in such samples under a magnetic field of 10T is 
about 10 -4 — 10~ 5 per nucleus, insufficient to gener- 



ate signal from a sample with 10 12 — 10 15 nuclei. [4] 
Higher levels of nuclear polarization are required for 
a detectible NMR signal. An alternative to standard 
NMR spectroscopy is optically pumped NMR spec- 
troscopy. In this approach non-equilibrium electron spin 
polarizations are generated through exciting electronic 
interband transitions. These spin-polarized electrons 
then polarize the nuclei during their relaxation. [4, 8-10] 
This process is known as dynamical nuclear polarization 
(DNP).[11] DNP has recently been achieved electrically 
in ferromagnetic-semiconductor heterostructures.[12] 

The nuclear spin coherence and relaxation dynamics 
will also be influenced by the fact that all nuclei in GaAs 
carry spins, so nuclear diffusion cannot be ignored. [9, 13] 
In the DNP process the initial spatial distribution of the 
nuclear polarization will follow the nuclear-electron mu- 
tual spin-flip rate, proportional to the square of the elec- 
tronic local density of states, for the hyperfine interac- 
tion is a point contact interaction. In lightly doped bulk 
semiconductors the initial polarization of the nuclei is re- 
alized around impurities, whereas in PQW's the confining 
potential creates a position-dependent electron distribu- 
tion, which will determine also the initial polarization of 
the nuclei through the sample. Under continual optical 
pumping the diffusion process will lead, eventually, to the 
uniform polarization of the nuclei in the whole sample. 
The value of the diffusion constant in GaAs systems was 
estimated previously by Paget [13] for bulk samples and 
by Malinowski and Harley[9] for square quantum well 
structures to be of the order of 10~ 14 — 10~ 13 cm 2 /s. 

Once polarized, the nuclear spin dynamics can be con- 
trolled directly with magnetic fields, or indirectly with 
electric fields. A method for electric field control of the 
nuclear spin relaxation time in quantum well (QW) struc- 
tures was first reported by Smet et al. [5] , and is based on 
electrically tuning the electron density in a QW across a 
Quantum Hall ferromagnet transition. In the low tem- 
perature regime such control is possible due to electric 
field induced modifications to the spectrum of collective 
excitations. 

In this work we analyze two different situations in 
which the nuclear spin relaxation time can be con- 
trolled with an external electric field. We treat the two- 
dimensional electron gas at much higher temperatures, 
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outside of the Quantum Hall regime. The DNP process 
remains very efficient up to 30K, and we predict a sen- 
sitive dependence of the nuclear spin relaxation time on 
the electric field, originating from the manipulation of the 
electronic local density of states across the PQW. This 
manipulation of the electronic local density of states pro- 
duces a position-dependent nuclear spin relaxation rate. 
The position dependence can be observed in a sample 
with a single 5-doped layer of an NMR active material 
inserted at a certain position inside the PQW. The in- 
serted layer needs to be made of different nuclei than 
the one present in the host PQW to ensure that its ob- 
served resonant frequency is not the same as that of the 
host PQW nuclei. Even in a sample without a specific 
(5-doped layer substantial control of the nuclear relax- 
ation rate is possible. For uniformly distributed nuclei in 
the PQW, the appropriately averaged nuclear relaxation 
time depends on the occupancy factor of the electronic 
conduction subbands, decreasing stepwise as the electron 
density increases. 

We report both analytical and numerical results for 
the electric field dependent nuclear spin relaxation time. 
The special structure of the confining potential in the 
PQW allows exact analytical solutions for the position 
dependent electronic density across the structure. More 
accurate results can be obtained from numerical calcula- 
tions of the electronic local density of states in the PQW 
based on a kp code, however we find that the analytic 
results are respectably accurate for the structures we con- 
sidered. We find that the electron and nuclear spin relax- 
ation times can be tuned over a range of at least an order 
of magnitude using an external electric field. To identify 
the most suitable PQW structure for the nuclear spin re- 
laxation time manipulation, we investigate two different 
possible samples, with different confining potentials. 

The paper is organized as follows. The second section 
is dedicated to the general presentation of the calculation 
of the nuclear relaxation time as a result of the hypcrfine 
interaction. The third section presents our analytical and 
numerical analysis of the PQW. Finally we summarize 
our results in the last section. 



II. THE GENERAL FORMALISM FOR THE 
ELECTRONIC AND NUCLEAR RELAXATION 
TIMES 

Our analysis of the electronic and nuclear spin relax- 
ation times in low dimensional systems follows the pre- 
vious calculation by Overhauser[14] for the case of three 
dimensional (3D) bulk metals. We only consider effects 
due to the hyperfine interaction which couple the elec- 
tronic and nuclear spins. The interaction Hamiltonian, 
known also as the Fermi contact term, can be written as 

H = # ( r «) =Ey (°™ • *e) 5{r - r «) > 

n n 

(2.1) 



where f} n and (3 e are the nuclear and electron spin mag- 
netic moments, a n and a e the Pauli spin operators for the 
nucleus and electron, and S(x) is the usual delta function. 
The argument of the delta function in Eq. (2.1), r — r„, 
represents the relative distance between the nuclear and 
electronic spins, and shows the local character of the hy- 
pcrfine interaction. The effects of the Fermi contact term 
are better understood if we express the product of the two 
Pauli spin operators, (a n ■ a e ), in terms of the creation 
and annihilation operators, c„( e ), as 

a n ■ a e = a z n a z e + 2 (a+a; + <j-<j+) , (2.2) 

an expression which clearly identifies the spin-flip pro- 
cesses involving both the electron and nuclear spin. The 
second term in Eq. (2.2) flips a nuclear spin from down 
to up in a process which implies also an electron spin flip 
from up to down. The last term in Eq. (2.2) describes 
the reverse process. 

To understand the spin relaxation process one has to 
consider the effects of a magnetic field 7i acting on the 
particle's spin. The final result of the interaction between 
the spin and the magnetic field will be an orientation of 
the spin, parallel or antiparallel, relative to the applied 
field. However, such a process does not occur instanta- 
neously, and its time dependence defines the character- 
istic spin relaxation time. The dependence on different 
parameters of the relaxation time originates from the par- 
ticular mechanism which produces the spin alignment. In 
our case this mechanism will be the hyperfine interaction. 
As the Fermi contact term involves both the electron and 
the nuclear spins we expect to obtain a single equation 
which will determine both the electron and nuclear spin 
relaxation times. 

The approximations we introduce are natural ones for 
nanostructures. We consider the electronic system to be 
in spin equilibrium with itself. Thus even though the 
electronic local density of states may vary from position 
to position, the electronic spin polarization is the same 
everywhere. Situations where the electronic spin system 
has position-dependent spin polarization can be consid- 
ered by combining the types of expressions found here 
with electron spin diffusion equations. For our purposes 
here we assume that at any specific time we have N + of 
the electrons with their spins oriented along the applied 
field (spin up) and N- with their spins oriented antipar- 
allel to the applied field (spin down). 

In contrast to the electronic system we permit the nu- 
clear system to develop spatially-dependent polarization. 
For simplicity we will derive the expressions for electronic 
and nuclear spin relaxation times considering only nuclei 
with total spin 1/2; the generalization for other nuclear 
spin values is straightforward. Thus for a given nucleus 
the probability of being spin up is M + (r„) and of spin 
down is M_(r„). One might be concerned if these local 
polarizations are well defined - it is known that in some 
circumstances non-Markovian effects can play a role, and 
we ignore such effects here. This approximation can be 
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interpreted in some circumstances as averaging over suffi- 
ciently long times, and in other circumstances of coarse- 
graining the nuclear system over a volume large com- 
pared to the internuclear distance, but small compared 
to relevant electronic length scales. This coarse-graining 
permits us to replace the orientation of a given nucleus 
in these equations with a statistical average. 

The time dependence of the relaxation process will be 
related to the time evolution of both the total electronic 
magnetization and the position-dependent nuclear mag- 
netization, which can be expressed, in terms of the differ- 
ences, D = N + - N_, and A(r„) = M+(r„) - M_(r n ). 

If we denote by (r„) the total numbers of electron 

spins which flip from down to up per second, and by 
W-+(r n ) the total number of electron spins which flip 
from up to down per second, each from interaction with 
nucleus n, we can write the change in D as 



dD 
~~dt 

and in A(r„) as 
dA(r ri 



dt 



2£[W + _(r„)-W_ + (r„)], (2.3) 



2[^_ + (r„)-VF+_(r„)] (2.4) 



The evaluation of (r„) and W ^(r n ) wu l be done 

based on Fermi's golden rule by treating the hyperfine 
interaction as a time dependent perturbation. 

Consider an electron transition between two states Iq 
and k'| in which the electron spin is flipped from down 
to up. According to Fermi's golden rule, the electronic 
transition probability from state k to state k' induced by 
nucleus n is given by 



2tt 



\H kk ,(r n )\ 2 N e (k\) 6(Ei - E f ) , (2.5) 



where -ffkk'( r «) = (^/(k' T ) \H(r n )\ ipi(ki)J represents 
the matrix elements of the hyperfine interaction Hamil- 
tonian for nucleus n, N e (k'^) is the electronic density of 
states for the final state, and Ei and Ef are the energies 
corresponding to the initial and final state. Accordingly, 
the total number of spin flips per second induced by nu- 
cleus n, W-\ (r„), is 



W + . 

E 

k 



W&(r n )M + (r n )/ FD (ki, Ek) [l - f FD (k' v E f )] , 

(2.6) 



where fpD(p<r,E) is the usual Fermi-Dirac distribution 
function for an electron with momentum p and spin a. 
The energy conservation expression 



2m 



— f3 e H — l3 n TL 



2m 



f} e U + (3 n H 



(2.7) 



allows us to eliminate kl and estimate (r„) using 

Eq. (2.6) by replacing the sum over the momenta with 



an integration over the energy using the density of states. 
The integration over energies can be performed simply 
in the case of small magnetic fields, f3 e ( n )TL <C fcsT (T is 
the temperature), for which the Fermi-Dirac function can 
be expanded in a Taylor series. W- + (r n ) is calculated 
in a similar way using M_(r n ) instead of M + (r„) and 
changing H to —H. 

After some algebra we find that the time dependence of 
the electron magnetization is described by the following 
equation 



dD _^ 512tt 3 /3 2 /3 2 N e \^(v n ) | 4 



~dt = ^ 



9hV 



(Do - D) 



^ 1024^/3^ e 2 A> e 2 fc B T|^(r ra )| 4 A 
■ 2^ (^o(r n ) - A(r„)) 



9h 



(2.8) 



where N e — J dsN e (e)f' FD (s), D and A (r„) are the 
thermal equilibrium values for D and A(r„), respectively, 
and ip(r n ) is the electronic wave function at the n'th nu- 
clear position. The first term on the right hand side 
(rhs) of Eq. (2.8) determines the electron spin relaxation 
time, whereas the second one the position-dependent nu- 
clear spin relaxation time. The generalization to the case 
of arbitrary nuclear spin / is straightforward if we con- 
sider that the nuclear spin flips probabilities W~ + (r n ) 
and (r„) describe transitions between adjacent nu- 
clear spin levels. 

In view of Eq. (2.8) the definitions of the electronic and 
nuclear spin relaxation times should be considered with 
some caution, as an exponential decay for the electronic 
magnetization is obtained only for the case when the nu- 
clear population approaches equilibrium, or is kept at a 
fixed noncquilibrium value. For the latter case a new ki- 
netic equilibrium value for the electron polarization, Dq, 
will be reached, which is determined by the value for D 
that makes dD/dt = in Eq. (2.8) for the constant 
noncquilibrium nuclear polarization A(r„). For calculat- 
ing the nuclear spin relaxation time we use the total spin 
angular momentum conservation 



dD_ 

dt 



2f(J+l)(2/+l) v dA(r„) 
3 2^' 



dt 



(2.9) 



assume the electron spin polarization is refreshed, and 
so remains approximately a constant, and separate the 
resulting equation into n equations: one for each nucleus. 

For the case of a two dimensional (2D) system, with 
a constant density of states, N e (e) = m e /2n (m e is 
the electron mass), the electronic and nuclear relaxation 
times arc 



1 v 512tt 2 m e 0j 0j |V>(r„)| 4 



and 



9h* (21 + 1) 



128tt m e /3 e 2 (i 2 n |V>(r»)| 4 fc 5 r 
3?i 5 7(/ + l)(27+l) 



(2.10) 



(2.11) 
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where A represents the area of the 2D system. Tem- 
perature corrections to these times are of the order of 
UbT/Ef, Ep being the Fermi energy of the electronic 
system. As we can see from Eq. (2.10) the electronic 
spin relaxation time is temperature independent (up to 
the first order of approximation) suggesting that the re- 
laxation mechanism due to the hyperfine interaction can 
be the dominant one at low temperatures, as the relax- 
ation times corresponding to other mechanisms increase 
as the temperature decreases. [14] According to Eq. (2.11) 
the nuclear spin relaxation time follows the well known 
Korringa law, which states that the ratio l/(T\ n T) is 
temperature independent. 

Let us turn our attention to the quasi-two-dimensional 
case of a QW. In such systems the growth-direction 
electronic motion is restricted to discrete energy lev- 
els through confinement in a potential well. Accord- 
ingly, the electron energy dispersion can be written as 
e(k) = Si + k 2 /2m e , where £j is the minimum of the 
conduction subband i and k is the electron momentum 
in the plane of the quantum well. If the separation be- 
tween the energy levels is large enough compared to the 
Fermi energy, the electrons will be frozen in the first en- 
ergy subband leading to the practical realization of a 2D 
electronic system. In this case the electronic and nu- 
clear spin relaxation times are given by Eqs. (2.10) and 
(2.11). However, if the Fermi energy exceeds the sep- 
aration between the energy levels, we must consider a 
multiply-occupied subband system. For the case of mul- 
tiple subband occupancy the density of states is 



N(e) 



' m e /27r, si < E F < s 2 
2m e /27r, e 2 < E F < e 3 

. lm e /2ir, ei < Ep < ei +1 



(2.12) 



where we consider that the l-th is the last subband of the 
system occupied by electrons. Two additional complica- 
tions arise that prevent direct use of this larger density of 
states in Eqs. (2.10-2.11). First, the envelope functions 
for the different subbands will differ, and second, elec- 
tron spin-flip scattering can take place inside the same 
subband or between two different subbands. 

The electron and nuclear spin relaxation times for a 
QW can still be written very simply, and physically, in 
terms of the electronic local density of states, 



T ^ ~ v S 



10247r 3 /3 2 /? 2 / deAj(r n ,e)f FD (e) 
V ^ m(2I+l)JdrdsA e {r,s)f' FD {e) 



(2.13) 



and 



T -i, x = 512^mkBTjdeAl(r n ,e)f FD (s) 
ln[n> 3W(J + 1)(2/+1) 

(2.14) 

The electronic local density of states at the nuclear posi- 
tion r„ is 

A e (r n ,e) = \^ m (r n )\ 2 5(£ - E m ) , (2.15) 



where m labels the state, E m its energy, and V , m( r «) its 
wavefunction at the n'th nucleus. 

In general for any QW structure the dispersion relation 
is quasi-two-dimensional, therefore, the electronic wave 
function will be written as a product between an envelope 
function, <j){z), and a Bloch function, u n K_(r), such that 
ipK,n{ r n) = cxp [iK ■ R]<fi(z)u nK (r n ), and 

A e (v n ,e) = ]T |^- K (r„)| 2 <5(£ - E jK ) , (2.16) 

where j is the subband index and K the momentum. 
The envelope functions for parabolic quantum wells will 
be evaluated both analytically and numerically. For the 
numerical estimations we will use a k • p calculation. [15] 
For the QW's stepwise density of states, 



/ 



deA 2 e (r n ,e)f FD (e) 



= | U (r„)| 4 V|^(r„)| 2 |0 fc (r„)| 2 9( £ , 



c{j,fc}) 



(2.17) 



with 0(sj) a temperature dependent factor. Such factors 
are important when the Fermi energy approaches an 
i.e., when the electrons start to occupy higher conduc- 
tion subbands. The energy integration over the step-like 
density of states in Eq. (2.17) leads to temperature de- 
pendent factors of the following form: 

= rr^-i • (2- 18 ) 



exp 



- E f 



k B T 



+ 1 



Eq. (2.18) clarifies the role of the temperature when 
higher subbands start to be occupied, i.e., the Fermi en- 
ergy is of the order of Ei. 

For K ~ the value of the Bloch function for the con- 
duction band can be approximated as a constant. The 
approximate constant value of the Bloch function at the 
71 Ga nucleus site, |w(r„)| 2 — 5.2 x 10 25 cm~ 3 , was ex- 
tracted by comparing our numerically calculated nuclear 
relaxation time with the available experimental data for 
a GaAs/Al .iGa . 9 As QW.[16] 

III. THE PARABOLIC QUANTUM WELL 

In this section we will discuss the nuclear spin relax- 
ation time for the case of parabolic QW's. The special 
form of the confining potential in these structures allows 
for exact analytical results which permit an evaluation of 
the nuclear spin relaxation time. The analytical results 
will be compared with numerical data obtained from the 
k • p calculation. 

An initial nuclear polarization obtained by DNP will be 
inhomogeneous, and for short times will be proportional 
to T^{z), so for one occupied subband the initial nuclear 
magnetization 



m(z,t = 0) cx \(j)(z)\ 4 



(3.1) 
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Our calculation assumes that DNP is performed in the 
absence of an applied electric field and for single subband 
occupancy. In this case (/>(z) corresponds to electrons 
situated in the first subband of the conduction band. 
The value of the nuclear magnetization will be influenced 
mainly by two processes, spin relaxation and spin diffu- 
sion. We track the evolution of the position dependent 
nuclear polarization in order to obtain the time depen- 
dence of the total nuclear magnetization in the PQW. 
Considering both spin diffusion and spin relaxation, the 
total nuclear magnetization across the PQW can be eval- 
uated as the solution to the following equation: 



dm(z,t) d 2 m(z,t) m(z,t) 



dt 



D 



dz 2 



T ln (z) 



(3.2) 



where D represents the nuclear diffusion constant. The 
combined effects of nuclear spin diffusion and relaxation 
indicate that the longer-time dynamics of the magnetiza- 
tion is non-exponential in PQW's.[16] If we consider that 
initially the PQW is in a state where only the first con- 
duction band is filled with electrons, during the DNP pro- 
cedure nuclei situated close to the sample's center become 
strongly polarized. The initial polarization is smaller for 
nuclei farther away from the center of the PQW where 
the electronic density is smaller. Even farther from the 
PQW center, where the electronic wave function is essen- 
tially zero, the nuclei are initially unpolarized. However, 
nuclear spin diffusion can bring the polarization gener- 
ated in the center of the well out to the entire sample. 

The position dependent initial nuclear spin relaxation 
time can be extracted from Eq. (3.2) if the nuclear dif- 
fusion term is absent (D = 0). As a consequence the 
observed total initial nuclear relaxation time is given by 



1 = JdzTr n 1 (z)V n (z) 
Ti„ JdzV n (z) 



(3.3) 



where V n {z) is proportional to the initial position- 
dependence of the nuclear polarization. For the case of 
optical pumping, where the nuclear polarization is real- 
ized through the hypcrfine interaction, V n (z) cx |0i(z)| 4 . 
Our estimation of the nuclear spin relaxation times will 
be done for temperatures where DNP in GaAs is efficient 
(T ~ 30K), however, we expect that the obtained results 
can be easily extended to higher temperatures. 



A. Analytical Results 

In an ideal parabolic QW (IPQW) the form of the 
confining potential is 



V{z) = ^m e uj 2 z 2 , 



where oj is the characteristic frequency. The electronic 
envelope functions and the characteristic subband ener- 
gies are obtained from the time-independent Schrodinger 



equation 



h 2 d 2 1 



2m e dz 2 2 



m e uj z | 0(z 



s) = e n (j)(z), (3.5) 



as 



<j> n {z) = iV„exp 



TO e W 2 
Z 

2h 



H n -i 



1/2 



and 



e n = [ n - - ) Uluq 



(3.6) 
(3.7) 



where the subband index n is a natural number, H n (x) 
arc the Hcrmitc polynomials, and N n is the normalization 
factor which can be obtained analytically if the width 
of the QW is much bigge r than the PQW length scale 
L 3> Iq (lo = \fhj (m e ojo)). According to Eq. (3.7) the 
energy distance between adjacent conduction subbands, 
AE = e n — e n -i, is a constant depending only on the 
confining potential's characteristic frequency u>o- 

An additional external electric field, J 7 , will change the 
Schrodinger equation by introducing a field dependent 
term of the form 



V(z) 



zTz 



(3.8) 



whose effects can be easily understood if absorbed into 
the parabolic confining potential. The external elec- 
tric field will not change the form of the envelope func- 
tions, however their positions in the QW are shifted 
by a distance proportional to the strength of the field, 
z = ef/m e WQ. The new wave functions are 



i(z,T) = A„exp 



(z - z f 



/m e u> \ 1/2 / \ 



. (3.9) 



The subband energy in the presence of an external elec- 
tric field will be lowered to a value 



e 2 T 2 
2m e Lu 2 



(3.10) 



where the energy difference, AE = e n — e n -\, is the 
same as in the absence of the electric field. The electric 
field will allow a direct control of the electron density 
profile across the sample leading to the manipulation of 
the nuclear spin relaxation times in PQW's. 

In the following we will consider two different PQW 
configurations, i.e, two structures with different charac- 
teristic frequencies lo : the first structure will be referred 
(3-4) to as IPQW I (u = 2.3xl0 13 s~\ AE = 15 mcV) and 
the second one as IPQW II (u = 7.6 xlO 13 s" 1 , AE = 50 
meV). This choice of the PQW's will allow us to discuss 
the relative influences of temperature, electric field, and 
well structure on the nuclear spin relaxation times. 



6 



0.03 



T 0.02 



C 



"0.01 



0.00 



— OkV/cm 


1 1: : 


- - 20 kv/cm 


-40 kV/cm 




< a > , ■ , A 


1 



0.10 



0.08 



« 0.06 
h 

= 0.04 
h 

- 0.02 



-600 -400 -200 200 400 

z(A) 



— kV/cm 
40 kV/cm 
80 kV/cm 



0.00 



-250 -150 -50 50 150 250 
Z(A) 




0.06. 



20 40 60 80 100 
F (kV/cm) 



FIG. 1: The position dependence of l/(Ti„T) for different val- 
ues of the external electric field [(a) IPQW I and (c) IPQW II]. 
The field dependence of the l/(Ti„T) ratio for 5-doped layers 
situated at different positions across the QW [(b) IPQW I 
and (d) IPQW II]. 



In Fig. 1(a) and 1(c) we present the position depen- 
dence of 1/ (Ti„T) as a function of the applied field for the 
case of single subband occupancy. IPQW I is more suit- 
able than IPQW II for electric field control, as the rate 
of change with field manifest in l/(T\ n T) is much larger. 
The displacement of the electronic wave function with an 
applied electric field is larger for wells with smaller values 
of the characteristic frequency, as z$ ~ uJq 2 , and IPQW 
I has a smaller u> than IPQW II. In Fig. 1(b) and 1(d) 
we plot the value of the 1 / (Ti„T) ratio as a function of 
the applied electric field for three different positions in 
the QW. These positions corresponded to initial nuclear 
polarization via DNP of the maximal value (center of the 
well) and half of the maximal value. The tunability of 
the nuclear spin relaxation times at these particular po- 
sition is a few orders of magnitude, with larger tunability 
for IPQW I. 

This property of the PQW's will allow direct control of 
the relaxation process in <5-doped layers of nuclei situated 
at certain positions inside the QW. The only requirement 
is that the <5-doped layer is situated in a region where 
the DNP procedure will produce a high enough initial 
polarization of the nuclei. The nuclei inside these layers 
will not be affected by nuclear diffusion, as this process 
requires similar nuclear species, or at least nuclear species 
with close values of the resonance frequencies. [11] 

In Fig. 2(a) and 2(c) we present l/(Ti„T) for uni- 
formly distributed nuclei as a function of the Fermi en- 
ergy (which depends on the electron density). Manipu- 
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FIG. 2: The initial value of l/(Ti n T) as a function of Fermi 
energy at several temperatures [(a) IPQW I and (c) IPQW II] . 
The electric field dependence of the initial nuclear spin relax- 
ation times for single subband occupancy, with and without 
nuclear spin diffusion [(b) IPQW I and (d) IPQW II]. 



lation of the electronic density in the QW can be accom- 
plished with a gate. The reference value for the energy is 
set at — AE below the minimum of the first conduction 
band. l/(T\ n T) depends stepwise on the Fermi energy 
at very low temperatures, increasing every time a new 
conduction electron subband becomes occupied. How- 
ever corrections to the Korringa law become important at 
higher temperatures for Fermi energies near these steps. 
As we can see from Fig. 2(a) at T=30 K for IPQW I, 
temperature-induced smearing will destroy the stepwise 
character of l/(7i„T). For IPQW II due to a larger ener- 
getic separation between subbands the stepwise behavior 
remains for higher temperatures. 

Fig. 2(b) and 2(d) present the initial value of the nu- 
clear spin relaxation time as a function of the applied 
electric field. The relaxation time, both in the absence 
of and presence of nuclear diffusion, is extracted as the 
first derivative of the total nuclear magnetization in the 
QW. For a nuclear spin diffusion constant close to the 
one extracted from observations in bulk GaAs systems, 
the diffusion process will diminish significantly the elec- 
tric field tunability of the nuclear spin relaxation time. 
Due to the larger distances the electronic wave function 
is displaced in IPQW I, the tuning calculated there is 
significantly more robust to nuclear spin diffusion than 
IPQW II. 

Thus we see that different PQW structures behave in 
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FIG. 3: The position dependence of l/(Ti„T) for different 
values of the external electric field [(a) PQW I and (c) PQW 
II]. The field dependence of l/(Ti n T) for 5-doped layers situ- 
ated at different position across the QW [(b) PQW I and (d) 
PQW II]. 



FIG. 4: The initial value of l/(Ti„T) as a function of Fermi 
energy at several temperatures [(a) PQW I and (c) PQW II] . 
The electric field dependence of the initial nuclear spin relax- 
ation times for single subband occupancy, with and without 
nuclear spin diffusion [(b) PQW I and (d) PQW II]. 



a complementary fashion in the presence of external elec- 
tric fields. If the electric field is used as an electron 
density gate, controlling the subband occupancy in the 
system, the stepwise dependence of 1/ (Ti n T) on electron 
density can be better observed in structures with a larger 
splitting of the energy levels, for larger subband separa- 
tions suppress the temperature dependence of the steps. 
On the other hand, if the electric field is used to tune the 
nuclear spin relaxation time by controlling the electronic 
local density of states (tilting the PQW), structures with 
a smaller value of the characteristic frequency, and im- 
plicitly with a smaller splitting of the energy subbands, 
are more resilient to nuclear diffusion. When the field is 
used to tilt the PQW, inserting a <5-doped layer of a dif- 
ferent material at a certain position across the QW will 
permit sensitive control of the relaxation time over a few 
orders of magnitude (see Fig 1 (b) and (d)). The <5-dopcd 
layer should be made of a different material, with a dif- 
ferent NMR resonance frequency than the host material. 



B. Numerical Results 

In this section we will numerically evaluate nuclear 
spin relaxation times for more realistic versions of the 
PQW's. The first structure (referred from now on as 
PQW I) is a 1000 A Al^Gai-^As system sandwiched be- 
tween two 100 A Alo.4Gao.6As barriers. [7] The value of 
the aluminum concentration across the QW varies from 



0.4 in the barriers to in the QW's center. The sec- 
ond structure (referred from now on as PQW II) is a 300 
A Al^Gai^ajAs system sandwiched between two 100 A 
Alo.3Gao.7As barriers. [6] In this case the aluminum con- 
centration across the sample varies from 0.3 in the bar- 
riers to in the center of the QW. Our estimation of the 
energy difference between the first and second conduction 
bands is AE = 16.6 mcV for PQW I and AE = 46.1 mcV 
for PQW II, values which are close to the ones consid- 
ered in the two analytic examples. The evaluation of the 
envelope functions in these particular structures is done 
based on a 14-band multilayer kp calculation. Therefore 
additional effects will be evident in the k • p calculations, 
including a diminished value of the energy difference be- 
tween higher conduction subbands. 

In Fig. 3(a) and 3(c) we present the position depen- 
dence of 1/ (Ti n T) for different values of the applied elec- 
tric field. As seen in the analytic studies, in PQW's the 
applied electric field leads to a displacement of the elec- 
tronic wave function with direct implications for the po- 
sition dependent nuclear spin relaxation time. Fig. 3(b) 
and 3(d) showcase the field dependence of l/{T ln T) for 
three different position across the QW's. These positions 
were selected in the same way as those for Fig. 1, corre- 
sponding to maximal and 50% of maximal initial nuclear 
polarization. The electric field tunability is much higher 
in PQW I, where the effects of the electric field are con- 
siderable even at low fields. The achieved tunability is 
a few orders of magnitude for both structures for any 
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of the three considered positions. The results shown in 
Fig. 3 are very similar to those obtained with the analytic 
model and shown in Fig. 1. In Fig. 4 we present the value 
of l/(T\ n T) for uniformly distributed nuclei as a func- 
tion of the Fermi energy. Just as for the analytic results 
presented in Fig. 2 the temperature-induced smearing is 
more significant for PQW I, which has a smaller subband 
separation, than for PQW II. However the effects of nu- 
clear diffusion are more pronounced for PQW I than for 
PQW II. 

IV. SUMMARY 

We have performed a detailed analysis, both analyti- 
cally and numerically, of the nuclear spin dynamics as a 
result of the hyperfine interaction in parabolic quantum 
wells. Analytical formulae for the electron and nuclear 
spin relaxation times were obtained for general 2D struc- 
tures and extended for quasi 2D structures. The elec- 
tronic structure enters these formulae in the form of the 
electronic local density of states. For the quasi 2D sit- 
uation we considered an ideal parabolic QW, for which 
analytic results are available for the field-dependent enve- 
lope functions. Thus we were able to explore the electric- 
field manipulation of the nuclear spin relaxation time 
entirely within an analytic model. These results were 
checked with numerical kp calculations for the envelope 
functions. The results obtained from analytical and nu- 
merical calculations are consistent, with small differences 
originating from approximations used in the analytical 
calculation. 

Our specific calculations were for two different PQW 
structures, where the main difference was the value of 
the characteristic frequency in the confining potential 



(or subband energy separation). Electrical field control 
of the nuclear spin dynamics in PQW structures can be 
achieved in two different ways: by controlling the elec- 
tronic subband occupancy or by controlling the position 
dependent electronic wave functions. When the elec- 
tronic Fermi energy is changed a stepwise dependence of 
the nuclear spin relaxation time is predicted. Our calcu- 
lation identifies the PQW with a larger energy splitting 
between subbands as more suitable for such manipula- 
tion. On the other hand, for single subband occupancy, a 
position dependent control of the nuclear spin relaxation 
time is possible as the electronic local density of states 
across the PQW is tunable by tilting the PQW potential 
with an electric field. In this case we suggested that the 
insertion of a <5-doped layer of different NMR active nu- 
clei than the host nuclei would allow manipulation of the 
nuclear spin relaxation time on the scale of a few orders 
of magnitude. Here samples with a small value of the en- 
ergy splitting between adjacent electronic subbands are 
more suitable. The effect of nuclear spin diffusion on the 
nuclear spin dynamics for single subband occupancy was 
considered for several values of the applied electric field. 
We find that the presence of spin diffusion will dimin- 
ish the electrical field control of the relaxation process. 
However, for the PQW with a smaller energy splitting of 
the subbands, electric field tunability of the nuclear spin 
relaxation time is more robust. 
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